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4ri ■ Abstract 

O ■ 

^_J \ The properties of the interface in a phase-separated solution of polymers with differ- 
ed ! 

^ ■ ent degrees of polymerization and Kuhn segment lengths are calculated. The starting 

■n^ ' point is the planar interface, the profile of which is calculated in the so-called 'blob 

a ■ 

O ' model', which incorporates the solvent in an implicit way. The next step is the study 

O ■ 

'—'I of a metastable droplet phase formed by imposing a chemical potential different from 

^ ■ that at coexistence. The pressure difference across the curved interface, which corre- 

^ \ sponds to this higher chemical potential, is used to calculate the curvature properties 

^ ■ of the droplet. Interfacial tensions, Tolman lengths and rigidities are calculated and 

Tij- ! used for predictions for a realistic experimental case. The results suggest that inter- 

O' 

cn ' faces between phase-separated solutions of polymers exhibit, in general, a preferential 

L| ! curvature, which stabilizes droplets of low molecular mass polymer in a high molecular 

K> ■ mass macroscopic phase. 



1 Introduction 

When two chemically different liquid polymers are mixed, in the majority of cases 
phase separation takes place, because of the low gain in entropy after mixing, easily 
overruled by an adverse enthalpy of mixing. The presence of a common solvent 
usually suppresses the drive to phase separate, but only when the solvent is present 
at volume fractions of typically 90% or more. When the solvent is present at a lower 
concentration, phase separation takes place, giving rise to phase regions differing in 



polymer composition and total polymer concentration. The driving force for the 
phase separation is a repulsive, possibly entropic, force between the polymers, or a 
difference in the quality of the solvent for the two polymers [I] . Phase separation in 
a polymer mixture can also occur such that the polymers reside in one phase and the 
solvent in the other (complex coacervation). 

Interfaces between phases of coexisting, thermodynamically incompatible polymer 
solutions are referred to as solvent-solvent interfaces. An example is phase-separated 
protein and (neutral) polysaccharide solutions found in food systems |2] . The solvent- 
solvent interface, in that case a water-water interface, is then situated between the 
coexistent protein-rich and polysaccharide-rich solutions. The interfacial tension of 
solvent- solvent interfaces is extremely low (typically a few /iN/m or less) so that the 
interface is highly deformable, e.g. by convective flows, and therefore difficult to 
investigate by classical methods (such as the Wilhelmy plate method). Methods to 
estimate the interfacial tension are capillary rise [3] , spinning drop [H El [71 E] , shape 
relaxation after deformation [HI |9l [10], or by manipulation with a laser beam |11] . 

The macroscopic phase separation takes place via a stage in which one phase is 
a dispersion of microscopic droplets in a continuous matrix of the other phase. This 
state might be called a solvent-solvent emulsion. Figure [U gives an example. Such 
an emulsion is often surprisingly stable. This may be because of the small density 
difference (both phases are typically 90% solvent), low interfacial tension or the exis- 
tence of a preferential curvature. In order to get a basic theoretical understanding of 
the properties that play a role in the stability of such solvent-solvent emulsions, we 
present calculations of the interfacial tension, preferential curvature and rigidity of 
equihbrium and quasi-equilibrated curved interfaces between polymer phases. This is 
done using the so called 'blob model' [12] . The blob model is based on a Flory-Huggins 
mean-field approach but it is adapted to include excluded volume fiuctuations which 
are correlated over a blob size (correlation length) ^. Since solvent-solvent interfaces 
are wide, it contains only weak gradients so that it seems appropriate to describe the 
interfacial region in terms of the van der Waals squared-gradient expansion. 

In this article, we discuss two simplified cases: (i) identical degrees of polymer- 
ization and solvent qualities, but allowing for a solvent gradient in the interface, (ii) 
different degrees of polymerization and solvent qualities (i.e. different Kuhn lengths), 
but neglecting the solvent gradient in the interface. Case (i) was previously treated 
by Broseta et al. [131 El US] and will be used as an introduction to case (n). A 




Figure 1: Water- water emulsion, from a phase separated aqueous solution of 10% fish 
gelatin (non-gelling), and 10% dextran (200 kDa). Light phase: gelatin-rich; dark 
phase: dextran-rich. 

preferential curvature is expected in case (ii) because the energy in the compositional 
gradient depends on the blob size ratio and the ratio of the number of blobs per chain. 

2 The blob model 



2.1 Approach by Broseta et al. 



The derivation of Broseta et al. [131 HH 113 lEj starts with the Flory Huggins model 
[17] for an isotropic polymer melt consisting of two polymer types of equal length A^. 
We denote the volume fractions of type A and type B as (J)a = 4> and 0^ = 1 — 0. 
Then, 
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where a is the dimension of the lattice which equals the Kuhn segment length. The 
Flory parameter x represents the interaction between the two polymer types. 

The Flory-Huggins model is used by Broseta et al. [13] for a mixture of two 
polymers well above the overlap concentration (semi-dilute solution) in a solvent 
which is a good solvent for both polymers. In that situation the blob model can 
be used to describe the polymer solution [12]. The blob model incorporates the 



effect of ffuctuations due to excluded volume interactions by replacing the Kuhn 
segment length a by a 'blob size' ^ and the polymer length A^ by the number of blobs, 
Nb = N/{cC,^) per polymer chain: 
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where K = 0.024 is a constant [13], [I5], related to the osmotic pressure, and where u 
represents the interaction between blobs of type A and B. The whole volume is now 
filled by 'blobs' of dimension ^ with each 'blob' either containing polymer of type A 
or type B. An essential feature in the blob model is that both the blob size ^ = ^(c) 
and the interaction parameter u = u{c) depend on the (total) monomer concentration 
c [121 [18]: 
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where the exponents z/= 0.60 and x= 0.22 are such that non-mean-field chain statis- 
tics is explicitly taken into account (x can be interpreted as the number of monomeric 
contacts between overlapping blobs). Furthermore, c* is the overlap concentration, 
Rg the polymer radius of gyration in dilute (non-overlapping) conditions, and Merit 
is the interaction at a reference concentration, for which we shall take the critical 
demixing concentration, Ccrit- 

The blob model for an isotropic polymer mixture can be extended to non-homogeneous 
polymer solutions. This is achieved by allowing (f) = (j){r) and c = c{r) to be position 
dependent and by adding squared gradient terms to the free energy: 
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This expression is used by Broseta et al. [14] to determine density profiles and surface 
tension of a planar interface. We revisit their analysis in the Appendix and point to 
some small errors in their original formulas. 

We can extend the above free energy to non- symmetric polymer solutions by 
allowing the blob size and number of blobs to differ for each polymer type. This may 
be due to a different degree of polymerization, different solvent qualities, different 



Kuhn lengths, etc. The difference in blob size {C,a and C,b) and number of blobs 
{Nh^A and A^6,b), then leads to the following modification of the expression for the 
free energy: 
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where we have defined an effective blob size as 
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The form of the free energy in Eq. (jH]) is derived by replacing ^ — )■ ^a/b and A^^ — )■ Nt,^A/B 
in the corresponding terms in the free energy in Eq.([S]). For the other terms in Eq.([5]), 
such as the term describing the interaction between blobs, the blob size is replaced by 
an effective blob size which should somehow be related to the two blob sizes C,a and 
C,B- We have chosen the expression in Eq.(I7]) to calculate the effective blob size since 
it reduces to ^ for the symmetric polymer system and because its value is dominated 
by the smallest of the two blob sizes. Naturally, other models for the free energy may 
be constructed and, in particular, one may propose a different expression to calculate 
the effective correlation length, but we expect that most of the physics involved is 
captured by the expression for the free energy in Eq.([6]). 

In the following sections, this free energy is applied to study spherically (and 
cylindrically) shaped droplets. In the next section, we assume that we can neglect 
the variation of the total monomer concentration profile in the interfacial region and 
assume that elf) =c everywhere. 

2.2 Calculation of the curved interface profile 

When we neglect the variation of the total monomer concentration, we have that 
c{r) = c and we may replace both the blob size and the interaction parameter by 
their constant values in the bulk region: ^ -^ C, and u -^ u (the bar denotes the 
corresponding bulk value). Furthermore, the free energy is a functional of 0(r) only 
and we can write for the grand free energy: 
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with the two asymmetry parameters a and tq defined as: 
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The parameter a represents the difference in degree of polymerization of the two 

polymers, while tq reflects the difference in Kuhn lengths, e.g. due to a difference 

in solvent quality or chain stiffness. The dimensionless interaction parameter u is 

defined by 
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which can also be written as: 
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where we have used Eqs.(jl]) and ( TTT1) and the fact that Wcrit = (1 + v^)^/(l + a). This 
last expression will be used to link our theoretical results to more experimentally 
accessible quantities. 

Next, we consider a spherically shaped droplet with (equimolar) radius R. The 
radius R is really a radius of curvature and its sign is chosen such that a positive 
curvature corresponds to the phase rich in polymer A (the 'liquid') residing inside 
the droplet with the (metastable) phase rich in polymer B (the 'vapor') outside. As 
a consequence, a > 1 means that the degree of polymerization inside the droplet is 
larger than outside, and ro > 1 means that inside the droplet the characteristic length 
scale (blob size) is larger than outside. The equimolar radius i? of a liquid drop is 
defined in the expression: 
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The excess (grand) free energy of the metastable critical nucleus is defined as: 
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with Pjj the vapour pressure outside the droplet and p£ = p^ + Ap is the hquid pressure 
inside (see remark below, however). A = 4:TtE? is the area of a droplet of radius R. 
The quantity cr(i?) is the surface tension of the droplet as a function of the radius R. 
This is the quantity that we wish to study and the above formula provides a way to 
determine it from Afi. 

In spherical geometry, the free energy in Eq.(|8]) is: 
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The Euler-Lagrange equation to minimize the above free energy is given by 
2m(0) 0"(r) = _l!!M ct)'{r) - m'(0)0'(r)2 + /'(0) - /i , 



where the prime denotes a derivative with respect to its argument. The procedure to 
determine (j{R) as a function of R is now as follows: 

(1) We first determine the bulk densities 0o/ and <^o,^ coexisting across a planar 
interface {R = oo), which we refer to as the the state of coexistence. The corresponding 
value of the chemical potential is then denoted as /Xcoex- Determination of /icoex, 
00,^ and 00,1; is achieved by solving the following set of equations for given ou and 
asymmetry parameters a and Tq: 

/ (00,1)) = A^coex , / (00,^) = /^coex , 

f{(f)0,v) - /icoex 0O,i> = /(0O,<?) " /^coex 00,^ = -Pcoex ■ (19) 

The resulting phase diagram is shown in Figure |2] for various values of a. 

(2) Next, we vary the chemical potential /x to a value off- coexistence. The consequence 
of choosing a value of /i different from ficoex is the coexistence of two pressures across 
a curved interface. In the simplest case of spherical droplets we have droplets of liquid 
(-R>0) for yU>yUcoex and droplets of vapour (i?<0) for yU<yUcoex- For given fi, u and 
a the densities 0£ and 0„ are determined from solving the following two equations: 

/'(0,)=/i, /'(0£)=/x, (20) 

with the corresponding bulk pressures calculated from 

/(0^) -/x0^ = -p„, f{(j)i)- iJ,(j)e = -pi. (21) 
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Figure 2: Polymer phase diagram of mixing for several asymmetry ratios a, repre- 
senting a difference in degree of polymerization (see Eq.l fTTI) ). w is a dimensionless 
interaction parameter given in Eq.( lT4l) . The solid circles denote the locations of the 
critical demixing points, c<;crit = (1 + V^)^/(l + '^i 0crit = l/(l + V^)- 

It should be noted that far outside the droplet (r — )■ oo) the density (or pressure) 
becomes equal to that of the bulk, 0(oo) =0i,, but that only for a very large droplet 
is the density (or pressure) inside the droplet (0(0)) equal to the bulk (0^). 

(3) As a next step, the Euler-Lagrange equation for 0(r) in Eq. (iT8l) is solved numeri- 
cally with the boundary condition 0(oo) = 0^. The resulting density profile 0(r) thus 
obtained is finally inserted into Eq. flTSj) to determine R and into Eq. flT7|) to determine 
!S£l and thus a{R). 

In Figure [31 two examples of o"(-R) are shown as a function of 1/R. These quantities 
are shown in reduced units which are such that all energies are in units of k'QT/Q{6u)2 
and all lengths are in units of Doo = ^/{Qu)2 , which is the interfacial width at infinite 
incompatibility of the two polymers, i.e. at infinite degree of polymerization |14j . 
This means that the surface tension and radius in reduced units are: 
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The determination of cr(i?) is quite elaborate and most of the time we are only inter- 
ested in its general shape. For this reason we perform a curvature expansion in 1/R 
in the next section to determine the parabolic approximation to cr(i?) shown as the 
solid line in Figure |3l 
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Figure 3: Surface tension of a spherical droplet as function of 1/-R, for ratios of degree 
of polymerization a = 1 and 8. For the reduced units, see Eg. (1221) . R>0 corresponds 
to droplets of the phase rich in polymer A and R<0 corresponds to droplets of the 
phase rich in polymer B. We have set a; = 4 and tq = 1 (equal blob sizes). The solid 
lines are parabolic approximations to o"(i?) determined from the curvature expansion. 

2.3 Curvature expansion 

We can analyze the results in the previous section by performing a curvature expan- 
sion in 1/R. All quantities are then expanded in 1/R. For the density profile and 
surface tension, for instance, we have: 

0(r) = 0oW + ^ + ... 
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25a 2k + k 
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where z = r — R and where a is the surface tension of the planar interface, 6 is the 
Tolman length yL9|, related to the radius of spontaneous curvature as Rq = 2k /5a 
[20] , k is the bending rigidity and k is the rigidity constant associated with Gaussian 
curvature [21]. An expansion of the Euler-Lagrange equation in Eq. (TT8|) yields to 
zeroth order: 

m(0o) 0o(^)^ = /(0o) - /^coex 0o(^) + Pcoex • (24) 

Expansion of the Euler-Lagrange equation to first order gives: 



2m(0o)0i(^) = -4m(0o)0o(^)-/^i-2m'(0o)0o(^)0'i(^) 

+<Pi{z) fr(0o) - 2m'(0o) 0o(^) - m"(0o) U4' 



(25) 



with yUi = 2cr/(0o,£ — 0o,«) [221 123]. One may show [22j that the coefficients in the 
expansion of a{R) in Eq. fl23|) can be written in terms of the density profiles (j)f){z) and 
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The surface tension and Tolman length are independent of the choice for the location 
of the interface (the position of the 2 = plane) but k and k do depend on this choice. 
For all of the above expressions we have taken the egumo/ar radius for i?, see Eq. flT5|) . 
which implies that the planar density profile obeys 



dz z (I)'q{z) = Q . (27) 

-J 

The same expansion can be carried out for a cylindrical interface. For the density 
profile and surface tension one then has instead of Eg. (1231) : 

0(r) = Mz)+^^+■■■, 

, , 5a k , , 

This gives us the opportunity to determine k and k separately with the result [22] 

GO 

k = fdz 2m{(j)o) z^ (p'oizf , 
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(29) 



The procedure to determine the curvature coefficients a, 6, k and k is now as follows: 

(1) The planar profile (f)o{z) is first determined from the differential equation in 
Eq. fl24l) with 0o,^, (po,v, /^coex andpcoex determined from the set of equations in Eq. flT9|) . 

(2) The location of the z = plane is chosen such that Eq. (l27|) is satisfied. This 
location depends on the asymmetry of the interface profile. 

(3) With 00 (-2) determined, a, 6 and k can all be evaluated from the integrals in 
Eqs.(l26D and (129]). 

(4) The constant fii is subsequently determined from fii = 2a/{(f)Q^i — (f)Q^y) which 
allows us to determine the bulk density values (pi//v from (pi//v = fJ'i/ f"{4>o//v)- 
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Figure 4: Variation of the Tolman length (in nm) as a function of monomer concen- 
tration for two ratios of the degree of polymerization and three ratios of the blob 
size. 

(5) For given 0o(-2) and /ii, the differential equation for 0i(2;) in Eg. (125!) may now be 
solved with the boundary conditions 0i(— oo) = 0£ and 0i(oo) = (j)v One may verify 
that if (j)i{z) is particular solution of Eg. (1251) then also 0i(-2) + constant ■ 0o(^) is also 
a solution. It turns out that k is independent of the value of this integration constant 
so we may arbitrarily take 0i (0) = for convenience [23] . 

(6) Finally, with 0i(2;) determined, k (or 2k + k) can be evaluated from the integral 
in Eg.(|29|). 

The result of this procedure for the Tolman length is shown in Figure H] for various 
values of a and tq. For a perfectly symmetric profile {a = l and ro = l) we have that 
(5 = identically: there is no preferred curvature to either phase. When a > 1, we 
obtain negative values for 5 indicating that the droplets of the phase rich in A (i?>0) 
have a higher surface tension as compared to droplets rich in polymer B (i? < 0). 
This feature can also be seen from the graph in Figure [3l The influence of tq is less 
pronounced with tq > 1 leading to a positive contribution and tq <1 leading to a 
negative contribution to 5. 

In Figure [5l we show the rigidity constants k and k and the combination 2k + k 
in units of k^T for a = l and for a = 8. It is noted that 2k + k is negative indicating 
that cr{R) as a function of 1/R has a negative second derivative, cf. Figure [3l 
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Figure 5: Variation of the rigidities (in units of kBT) as a function of monomer 
concentration for two ratios of the degree of polymerization, k is the bending rigidity 
and k is the rigidity constant associated with Gaussian curvature. 

3 Discussion 

In this section, we discuss in more detail the symmetric polymer system for which the 
actual numbers quoted in Figures H] and [5] are derived. The bar to denote bulk values 
shall be omitted in this section. We first consider the symmetric, planar interface as 
treated by Broseta et al. [TSj [HI [T5] . 



3.1 Symmetric polymer system 

For the symmetric polymer system, the influence of the solvent on the interfacial 
tension of the planar interface can be investigated using the approach of Broseta et 
al. [131 HH HH], detailed in the Appendix. In order to make the connection with 
experiments, we consider the parameters typical for aqueous solutions of gelatin and 
dextran [16] , which have a radius of gyration R^ ~ 18 nm, degree of polymerisation 
A^ ~ 1000, and molecular mass of the monomers, Mmon — 120 g mol~^. The overlap 
concentration c* may then be estimated from 

which gives c* ~ 0.82% mass per mass solution. By convention, we shall write con- 
centrations in units of % mass per mass solution which is achieved by multiplying the 
concentration expressed as number of molecules per volume by the molecular mass 
and dividing by Avogadro's number and the solvent (water) mass density. 

The blob size ^ in this symmetric polymer system can be calculated as a function 
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of the bulk monomer concentration c using Eq. (|3]) [18] : 

e(c)^0.43i?g(^-j . (31) 

From the blob size, the number of blobs per chain as a function of the total monomer 
concentration is then determined from: 

^b(c) = — ^. (32) 

For the polymer solution consisting of gelatin and dextran, the critical demixing 
concentration is estimated as Ccrit — 3.5% mass per mass solution. The blob size at 
the critical concentration can then be estimated from Eq.f pT]) . giving ^^rit — 2.6 nm, 
and the number of blobs per chain at the critical demixing concentration is estimated 
from Eq.(|32D as A/'b,crit ^ 325. 

By choosing in Eq.(|l]) the critical concentration of demixing as the reference con- 
centration, u also becomes experimentally accessible: 

^(^) = i;j— [— ' (33) 

where we have used that uj = N\yU [T^. For the symmetrical polymer system (a = l), 
we have that Wcrit = 2, which gives Merit — 0.0061. In Table 1, we list the values 
of the parameters of the gelatin-dextran system together with the parameters calcu- 
lated from these numbers and their range when the monomer concentration is varied 
between c= 3.5% (w/w) and c= 10% (w/w). 

The free energy by Broseta et al. in Eq.(|5]) for the symmetric polymer system 
takes into account the reduction of the monomer concentration by the accumulation 
of solvent at the interface, which is ignored in the calculation of the properties of non- 
symmetric interfaces. It is therefore useful to first consider the relevance of this effect. 
This can be done by plotting the ratio of the interfacial tension of the planar interface 
calculated with and without taking the reduction due to the solvent redistribution 
into account. This ratio is (see Eq. fjiOj) ): 

Q ^ (yC=C{z) ^ 1- A1-MA2 _ ^g^^ 

^c=constant ^ ^1 

It turns out that the effect of solvent redistribution is significant in the whole con- 
centration range, reducing the interfacial tension by at least 20% (See Figure [6]). 
However, the reduction is nearly constant above a monomer concentration of 1.5 Ccrit 
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quantity 


value 


Rg 

N 


radius of gyration 
degree of polymerisation 
monomer molecular mass 


18 nm 

1000 

120 g mol"^ 


C* 


monomer overlap concentration 


0.82% (w/w) 


Ccrit 


critical demixing concentration 


3.5% (w/w) 


Scrit 

^b,crit 

Wcrit 


critical blob size 

critical number of blobs per chain 

critical interaction parameter 


2.6 nm 

325 
0.0061 


C 


monomer concentration 


3.5-10% (w/w) 


iVb 


blob size 

number of blobs per chain 

interaction parameter 


2.6-1.2 nm 

325-1200 

0.0061-0.0084 



Table 1: Values of the parameters of the gelatin-dextran system [16] and their varia- 
tion with monomer concentration c. 




Figure 6: Ratio of the interfacial tension corrected and not corrected for the effect 
of solvent redistribution as a function of the (bulk) monomer concentration. In this 
example, we have set u= 0.05. 
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{uj > 3.7), above which the decreasing amount of solvent is compensated by the in- 
creasing segregation of the polymers at the interface. It is tentatively assumed that 
calculations which ignore the solvent become more reliable at increasing monomer 
concentration. The properties of the curved interface should therefore be considered 
only at w > 4. 

3.2 Asymmetric polymer system 

In our theoretical analysis, the Tolman length and rigidity constants are calculated, 
as a function of the dimensionless interaction parameter u, in reduced units similar 
to those in Eq.l^: 

^^A = Mi5 and ~k^'-Plk. (35) 

To transform these units into experimentally relevant values, we need to make as- 
sumptions on the value of the effective blob size ^ (Eq.([7])) and the blob interaction 
parameter u (Eq. flT3|) ) also for the asymmetric polymer system. [The value and con- 
centration dependence of u is given by Eq.( !T4l) .] Since the concentration dependence 
of ^ and u is the same as in Eqs.( l3T]) and ( l33l) . we therefore only need to assume 
values for their critical values ^crit and Merit- As a first guess, it seems appropriate to 
assign for ^crit and Merit, the same values as in the symmetric case, ^crit— 2.6 nm and 
Merit ^ 0.0061. 

The Tolman length (Figure Hj) is, as obtained before [23], in the range of nanome- 
ters. The Tolman length is related to a preferential curvature of the interface and it 
can differ from zero because of unequal degrees of polymerization or because of un- 
equal blob sizes. For a prediction of the sign of the Tolman length for a specific pair 
of polymers, one therefore needs to know the molecular mass ratio and the respective 
Kuhn lengths. 

The implications of the sign of the Tolman length are apparent in Figure |3] which 
shows the interfacial tension as a function of radius of curvature for a symmetric and 
asymmetric polymer system. According to the definition of the asymmetry param- 
eter a in Eq.f jTTj) . the droplet phase has a higher degree of polymerization than the 
surrounding infinite phase when a> 1. Therefore, the shape of cr{R) as a function 
of 1/ R in Figure [3] implies that in a phase-separated polymer mixture, a droplet 
of high molecular mass has a higher surface energy than a droplet of low molecular 
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1/R 

Figure 7: Excess free energy density relative to the planar state AQ (in units of 
/cb^/(10 nm)^) as a function of the inverse radius (in units of nm^^). The solid 
line corresponds to droplets containing low molecular mass polymer surrounded by a 
continuous phase consisting of the high molecular mass polymers. The dashed line is 
the reversed situation. The maximum (solid circle) in the free energy corresponds to 
the preferential radius determined from Eq. fl38|) . In this example a = 8 and c/ccnt = 
1.8. 

mass. This would imply that the formation of droplets of the lower molecular mass 
is energetically more favourable. 



3.3 Droplet stability 

To elucidate the role of the radius dependent surface tension on the stability of the 
droplets, we consider the free energy of a dispersion of A'^s spherical droplets with 
radius R: 



n{Ns,R) = Ns4:nR^ 



25a 2k + k 



a 



R 



i?2 



A-TT 

-ApN,—R\ 



(36) 



The above free energy should be minimized with respect Ng and R keeping the total 
volume of the particles, Vo= Ng (47r/3)i?^, fixed. This yields the Laplace equation for 
the pressure difference: 



Ap 



2a 25a 



R i?2 



and we obtain for the preferential radius: 

R 



pref = 25+[4:5'--{2k + k) 
a 



(37) 



(3^ 



The role of the preferential radius is illustrated in Figure [71 Using the parameters 
listed in Table 1, the excess free energy density relative to the planar state, AQ = 
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Figure 8: Preferential radius (in nm) and corresponding excess free energy density 
relative to the planar state (in units of k^T j{\S] nm)'^). The solid line corresponds to 
droplets containing low molecular mass polymer surrounded by a continuous phase 
consisting of the high molecular mass polymers. The dashed line is the reversed 
situation. In this example a = 8. 

[n(i?) — n(i? = oo)] /Vo is plotted as a function of the inverse radius l/i?. The radius 
R is the droplet radius of the low molecular mass droplets (solid line) or the droplet 
radius of the high molecular mass droplets (dashed line). In either case, the maximum 
in the free energy corresponds to the 'preferential radius determined from Eq. fl38|) . 
The shape of the free energy shows the existence of metastable (and possibly long- 
lived) droplets with a preferential size that ultimately grow to form a single phase 
{R = oo). Figure[7|also shows that the system's free energy is significantly lower for the 
formation of droplets of low molecular mass polymer. This result indicates that when 
the emulsion is initially formed, the thermoydnamic path followed toward macroscopic 
phase separation is preferably through the formation of droplets containing polymer 
with the lower molecular mass. This effect should be observable by dynamic light 
scattering, when the two polymers have a high degree of monodispersity and the 
phases have a small difference in density. 

The preferential radius and the corresponding free energy are shown as a function 
of concentration in Figure [HI The preferential size of the (metastable) droplets is in 
the range of tenths of nanometers with the corresponding excess free energy of the 
polymer system of the order of a few k-^T per volume of (10 nm)^. 
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4 Summary 

The interfacial tension between coexisting incompatible polymer solutions is calcu- 
lated as a function of concentration for different degrees of polymerization and Kulin 
segment lengths. This is done using the blob model for semi-dilute polymer solutions. 
It turns out that, when the degree of polymerization or the Kuhn length of the in- 
compatible polymers differs, the interfacial tension is lowest when the interface has 
a certain preferential curvature. As a consequence, droplets of low molecular mass 
phase surrounded by a high molecular mass bulk phase may have an energy which is 
lower than that of the reversed situation. The preferential size was found to be of the 
order of tens of nanometers, and therefore observable by light scattering. 

A Original analysis by Broseta 

In this Appendix, we follow the analysis by Broseta et al. [131 [IH US] of the profile of 
a planar interface between two coexistent polymer phases with solvent. The polymers 
have equal degrees of polymerization and are chemically identical with respect to the 
solvent. We start with expression for the free energy in Eq.Q for a planar interface: 



dz 



3 '<^) + %:# 1^(1 



e 24^0 24^(1-^ 



^i^.a 



+ 



d{zf 



24^c2_ 

Broseta et al. [ll] continue by defining the dimensionless 77(2;) and e{z) as: 

_ l + ri{z) 



and 



ciz 



c + uceiz) . 



(39) 



(40) 



The free energy per unit area then becomes 



Ak^T 



dz 
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^3 24^(l-r/2) 24^(1 + Me)2 



(41) 



Next, an expansion is made by Broseta et al. ^14j assuming that e{z) is small every- 
where. One then has: 



N,i^ 



M' 



;i + ue) 
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(42) 



To second order in e, we then have for the free energy 
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2(3i/-l)2 

where we have defined u = NbU [1_4J and omitted constants and terms proportional to 
a constant times e. These are adsorbed in a constant bulk pressure p and chemical 
potentials /i^ and fie- As a final step, we introduce x = z/Doo, with Doo =^/{6u)2 
and the excess grand free energy due to the interface becomes: 
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If we compare the expression above to (A10)-(12) in Broseta et al. [H], we see that 
the prefactor of e'(x)2 should read u^ and not u. 

The above free energy is minimized in two steps. First, the profile ^0(2;) is deter- 
mined assuming that e{x) = 0. One finds that the bulk value f/ [ri^ = fi and rj^ = —ff) 
is determined by (/x^ = by symmetry) : 






^ = 0, 



(45) 



which has solutions for uj>uJc = 2 (77^ = 0). The profile riQ{x) is determined by solving: 
r]'^{xf (l+r^o) ,^,^ _^, , (l-r^o) 



2 
ln(l-r7o)-^+p. 



(46) 



4(1-77^) 2a; ""'" ' '"" ' 2uj 

Second, using rjQ^x), the profile £0(2^) is determined. In the bulk £ = which leads to 
the following expression for /i^: 
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:i + '')m+.,) + (Lz^i„(i_,) 
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The profile eo{x) is determined by solving: 
u 



2 4'(^) 



l+''«),„(l+,,„)+(l^l„(l-„„) 
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Again, comparing with (18) in Broseta et al. 114], we see that the additional factor 
of u in the expression for the free energy in Eq. (H^ leads to the presence of a factor 
u in front of Sq^x). 

With the help of these two profiles, the surface tension can then be calculated as: 
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(49) 



(50) 



One sees that Ai is straightforwardly calculated from the above integral, but that for 
the calculations of A2 one needs to determine both density profiles riQ{x) and e{x). 
The expression for A2 also differs from the one given in (29) in Broseta et al. 
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